Quaternions and hybrid nematic disclinations 
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Abstract 

Disclination lines in nematic liquid crystals can exist in different geometric conformations, char- 
acterised by their director profile. In certain confined, colloidal and even more prominently in chiral 
nematics, the director profile may vary along the disclination line. We construct a robust geomet- 
ric decomposition of director profile variations in closed disclination loops based on a quaternion 
description and use it to apply topological classification to linked loops with arbitrary variation 
of the profile. The description bridges the gap between the known abstract classification scheme 
derived from homotopy theory and the observable local features of disclinations. We compare the 
resulting decomposition of disclination loop features to a similar decomposition of nematic textures 
on closed surfaces. 



I. INTRODUCTION 



The defects in condensed materials have a significant effect on their physical properties. 
While their existence may be caused by impurities and flawed preparation process, controlled 
manipulation of defects can enable tuning of the physical characteristics of the material. The 
ability of a medium to host defects is a side-effect of symmetry breaking on the microscopic 
scale. As such, the defects are intimately related to the topological properties of degenerate 
microscopic degrees of freedom. In materials without broken translational symmetry, such 
as uniaxial and biaxial nematics pQ, ferromagnets [2] and superfluid helium [3111], the defects 
can be roughly classified by their dimension into point defects, line defects and walls [5]. 
Their classification with the aid of homotopy groups is a well-established and efficient way 
of describing their interactions and the conservation rules they obey [6]. However, the 
topological analysis does not account for the fine geometric details, which are under control 
of the free energy. 

Specifically for nematic braids - networks of closed nematic disclination lines [7HH] - recent 
development has shown that additional geometric constraints allow for a finer classification 
of nematic braids and simplifies the understanding of their rewiring [101 E] . The constraint 
of a fixed —1/2 disclination profile that holds for nematic braids is not universal and does 
apply necessarily to some constrained [12] or chiral systems [T3HT5] . But even without 
constraints, it should still be possible to convey the topological information through the 
geometric primitives. In this paper, we demonstrate a formalism that can describe the most 
general types of nematic disclinations in terms of the local director behavior and enumerate 
them using quaternions. We show that our interpretation naturally extends the existing 
formalism for —1/2 disclination loops and reveal a connection between closed disclination 
loops and two-dimensional nematic textures. 

II. OVERVIEW 

Before introducing any new results, we shall briefly overview the existing theoretical 
background regarding nematic defect loops. 

A nematic liquid crystal consists of elongated molecules that are locally orient at ionally 
aligned [I]. In a continuum approximation, the order is described by the director n, a unit 
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vector field, pointing in the direction of average molecular orientation. The head-to-tail 
symmetry of the molecules causes the sign of the director to be ambiguous. 

Topologically, the director field is a map from the coordinate space R 3 to the real projec- 
tive plane RP 2 - the ground state manifold (GSM). The fundamental group tti(RP 2 ) = Z2 
distinguishes disclination lines, where the director makes a half-turn around the disclination, 
from nondefect states, but makes no distinction between different disclinations [6]. 

Our specific interest is focused on disclination lines that are closed into a loop. Disclina- 
tion loops may be linked by other disclinations or carry a topological point charge, measured 
by the second homotopy group. Closed loops therefore also carry an additional topological 
index v e Z 4 . An established result shows [T6HT8] . that a set of n disclinations obeys a law 
in the form 



where q is the topological point charge of the set of disclination loops, Lkjj are linking 
numbers between loops and the indices i and j run over all the loops. Odd indices v 
correspond to disclinations that are threaded by other disclinations an odd number of times 
while even indices correspond to those that are threaded by and even number, which also 
applies to unthreaded loops. This theory describes all existing topological properties of 
disclination loops. However, the index v is an abstract attribute assigned to each of the 
loops, so the application to the systems that are found in experiments, simulations and 
theoretical models requires a way to relate this index to observable attributes. 

If instead of a three-dimensional space, the molecules of the nematic are restricted to 
move and orient in a two-dimensional plane, the system exhibits point defects, which can 
have any half-integer winding number, the lowest ones being —1/2 and +1/2 [6]. These 
two-dimensional point defects represent a restricted set of cross sections the disclination 
lines in three dimensions can have. The topological equivalence of all disclination lines in 
the three-dimensional case is caused by the ability of the director to rotate through the third 
dimension. If the director is forced by the free energy minimization or geometric constrants 
to be perpendicular to the disclination tangent, which is excellently obeyed in well-researched 
nematic colloids P [191 ED], closed —1/2 disclination loops can be described as ribbons that 
follow the rotation of the profile. Ribbons can be assigned an invariant called the self-linking 
number SI (211 E2], which in the case of —1/2 disclinations quantifies how many times the 
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three-fold disclination profile rotates around the tangent when we traverse the loop [10]. In 
a braid of many disclination loops, each self-linking number is of the form SI = m/3, where 
the numerator m is odd if the loop is threaded by an odd number of other disclinations, and 
even if the loop is linked by nothing or by an even number of disclinations. The self-linking 
numbers of a set of n disclination loops with a —1/2 profile obey a rule 



that bears resemblance to the Janich's rule ([T]). The integers Lkjj are the linking numbers 
between loops and Slj are self-linking numbers of individual loops. The comparison of 
equations ([!]) and (|2j) suggests a connection between the index v and the self-linking number. 
In the following section, we confirm and extend this notion to all disclination loops. 

III. THE PROFILE CIRCLES AND QUATERNIONS 

All disclinations in a 3D nematic have in common that the director on a test loop around 
the disclination traces an irreducible path in the ground state manifold (GSM). As the 
GSM is a sphere with additional condition of the antipodal points identified, the loop is 
topologically equivalent to a half of a great circle, which can be symmetrically extended to 
the full circle (Fig. [TJ. As the vicinity of the disclinations is subject to elastic free-energy 
restrictions, this is a quantitatively reasonable approximation. 

A closed disclination loop can be wrapped in a torus, S 1 g) S 1 , generated by revolving a 
disclination-encircling loop parameterized by u, around the disclination loop parameterized 
by v. For each value of v, the parameter u traces a different great circle in the GSM (Fig. Ua). 




Each circle is fully specified by its normal and a reference director at a chosen parameter, 
u = (Fig. [l]3,c). The sense of circulation around the circle is important and prescribes the 
sign of its normal. For example, the —1/2 profile corresponds to the normal of the circle 
being anti-parallel to the disclination tangent in the real space, as the director rotates in 
the opposite sense compared to the circulation of the enclosing loop (hence the minus sign 
in —1/2). The circle is invariant to rotations for 7r around its normal because each point in 
the real space maps to two antipodal points in the GSM. 

Let the parametrization of the director field on the torus be n(u, v). When v changes, the 
circle rotates, so the director at all values of parameter u transform with the same rotation, 
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(a) (b) ground state manifold 




relative rotation R'(v) changes along the loop 



FIG. 1: (a) The director profile specified on a torus enveloping a general disclination loop. The 
small radius is parametrized by u and the large radius is parametrized by v, producing a seamless 
doubly periodic parametrization. (b) The director on any small cross section lies on a great circle 
in the GSM. The origin of the parameter u = provides a reference in the plane of the circle 
and the sense of rotation in GSM with increasing u orients the circle's normal so that it defines 
a clockwise motion according to the right-hand rule, (c) Each circle in the GSM matches a cross 
section in real space with a fully specified parametrization u. (a-c) The orientation of the director 
cross section and its associated circle changes along the disclination with parameter v, which is 
described by a continuously changing rotation R'(v). 

parameterized by v, 

n(u,v) = R(v)n(«,0). (3) 

The n(u, 0) is simply a reference profile at a freely chosen point on the loop. This sets the 
initial rotation to unity, R(0) = 1. 

The rotations vary continuously with loop parameter v and thus describe a path in the 
space of all rotations in three-dimensions. We will use the quaternion representation of the 
SU(2) group. The unit quaternions 8, j and k obey the conventional structure relations 
d 2 = j 2 = k 2 = — 1, d| = —jo = k, jk = — kj = 5 and k§ = —5k = j. A rotation around an 
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arbitrary axis a by an angle follows the usual spinor formula 

R = cos ^ + (a x 5 + a y l + a z k) sin ^. (4) 

Rotation by 2tt around any axis yields R = — 1, which is the irreducible rotational motion 
that ends up in the same state where it started. Note that if we took the SO (3) representa- 
tion, these irreducible paths would be considered equal to the identity transformation, and 
some of the topological information would be lost. 

The continuity condition for a closed loop is satisfied if the director is continuous for all 
u, that is, if the circle we obtain after completing the loop at (v = 2tt) matches the circle 
at the beginning of the loop (v = 0) up to a it rotation around the loop tangent, which we 
align with the z-axis. Because the SU(2) group covers the space of rotations twice, the set 
of possible cumulative rotations expands to 

R(v = 2tt) G {l,k,-l,-k} = V. (5) 

The allowed rotations form a cyclic group Z 4 , enumerated by an integer u, which can be 
shown to correspond to the index introduced by Janich ([I]). As R(27r) is simply the total 
rotation of the director at any fixed u on a circuit around the disclination, it measures the 
fundamental group of the nematic encircled by the disclination. R = ±k are nontrivial 
rotations by 7r and mean the disclination is threaded by another disclination. R = 1 and 
R = — 1 correspond to an unlinked loop with even and odd topological charge, respectively 
(this topic is discussed thoroughly in [TS]). 

Each rotation can be decomposed into the rotation Rq{v) of the local coordinate frame 
relative to the global coordinate frame, and the rotation R'(t> ) of the director relative to the 
local coordinate frame. 

n(u,v) = R (v)R'(v)n(u,0). (6) 

We set the local coordinate frame to keep its z-axis aligned to the disclination loop tangent 
and vary continuously along the loop (e.g. the Frenet-Serret frame). We also ensure our 
framing is topologically equivalent to a framing of a planar circular loop without torsion, 
so the local coordinate frame simply rotates by 2tt when we traverse the loop: Ro(27r) = 
— 1. This extra rotational offset is present for all the loops, so we can factor it out. The 
multiplication of k u by —1 simply shifts v by 2 so an alternative parameter v* = v + 2 
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mod 4 can be introduced, rewriting the Janich's law Q as 

/ n n \ 

2 E^*- 2 E Lk y = 9 mod2. (7) 

V t=l i>j J 

We can see that the number of loops n plays a role in the conservation law because a trivial 
disclination loop with v* = holds a nonzero topological charge (one example is a Saturn 
ring defect). All the significant information is now encoded in R'(v), which is measured in a 
convenient disclination-aligned local coordinate frame. 

IV. ALGEBRA OF PROFILE TRANSITIONS 

Each geometric operation on the disclination profile can now be directly translated into 
the quaternion language. First, we reproduce the result for the —1/2 disclination profile. 
In the relative frame along the curve, the director lies in the xy plane and can only rotate 
around the z-axis, as it has to stay perpendicular to the disclination tangent. Rotation 
of the profile consists of both the passive rotation of the coordinate frame and the active 
rotation of the director. As the director lies in a plane, it can be specified by an angle a, so 
the profile rotation by ip yields a (f) = 3/2-0 rotation of the director, and subsequently the 
entire circle (Fig. |2^,b): 

11 3 

a (u) = --u; a'(u) = -^(w- VO +V> = at(u) + -ip . (8) 

In a loop with a self-linking number SI, the profile rotates by a total angle of ip = 2n SI, 
and the equation Q gives 

„„ N 3/2(2vrSl) n 3/2(2vrSl) „ SS] 

R'(2tt) = cos ^ + ksin^^ '- = k 3S1 . (9) 

Comparison to k u * yields the index v* = 3S1. The factor of 3 elegantly removes the frac- 
tional quantisation of the self-linking number that stems from the three-fold symmetry of 
the profile. This result immediately transforms the expression ([7]) to ^ and explains the 
correlation between the linking of disclinations and the self-linking number, discussed in 
Ref. pig. 

A similar calculation yields <fi = ip/2 and v* = SI for disclinations with a +1/2 profile. 
There is more to disclinations than just rotations around the tangent. An obvious pos- 
sibility is rotation around a perpendicular axis for <p = Pit, where P = ±1 is the parity 
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FIG. 2: (a) Notation used in equation ([8]). The —1/2 profile is two-dimensional, so the director 
can be represented by a single angle a that varies with parameter u. (b) Rotation of the director 
around the z axis induces a rigid rotation of the —1/2 disclination profile. Depicted is a rotation 
of the director for <f> = w, which yields a 1/3 of a full turn for the profile {ip = 2n/3). (c) The 
transition between +1/2 and —1/2 disclination profiles is realised by a rotation for it around a 
perpendicular axis (x axis, marked in black). There are two senses of rotation - the one depicted 
has the parity P = +1. 

that distinguishes left-handed and right-handed rotation for 7T. The result of such a rotation 
is the transition between +1/2 and —1/2 defect profiles (Fig. [2^). The axis of rotation in 
equation Q can point at any angle 6 in the xy plane, a = (cos#, sin#, 0), which leads to 

R' ± = P(8cos0 + jsinfl) = Pnexp(-[k0). (10) 

The rotations can be composed sequentially in the order they are encountered with in- 
creasing parameter v. A real disclination loop can have a profile that rotates around a 
gradually varying axis that is neither parallel nor perpendicular to the disclination tangent. 
However, the principal reasoning in our derivation is topological, so the exact geometric re- 
alisation of the rotations is not important. For example, the location of a transition between 
profiles can be pinpointed to the place where the circle's normal is exactly perpendicular 
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to the disclination's tangent. Whether the transition between — 1/2 and +1/2 disclination 
profiles happens along a very short segment of the disclination, or takes the entire disclina- 
tion length to complete, has no effect on the calculated index, as the director field in the 
second scenario can be smoothly combed to the first one without changing the topology. 
Consequently, rotations around the z axis (R' z ) and rotations for tt around a perpendicular 
axis (R'jJ are enough to describe all director profile variations: 

R' 2 = cos — + ksin — = exp(k|), (11) 
R' ± = P(ocos# + jsintf) = Poexp(-ktf). (12) 

For the rotations R' z , the angle is the local rotation angle of the director, which relates to 
the physical rotation angle ip of the profile: <fi = (3/2)-0 for a —1/2 profile and <fi = (l/2)ip 
for a +1/2 profile. 

Because of the presence of the unit quaternion 5 in the perpendicular rotation, it does 
not commute with z-axis rotations. 

R'MR' ± [e] = R'AOKi-^] (is) 

The rotation of +1/2 has the opposite effect than the rotation of —1/2 profile. The product 
of all local rotations must satisfy the continuity condition ([5]). As a perpendicular rotation 
introduces a factor 8, they must always appear in pairs - if the disclination turns from —1/2 
to +1/2 it must turn back into —1/2. The offset in the axis orientations adds an extra 
amount of torsion, 

#±[Pi,0i]# ± [P2,e 2 ] = -PiP 2 ex P (k(e 1 - 2 )). (14) 

If two sequential rotations about the same perpendicular axis are performed with the same 
parity (twice by tc or twice by —tt), they yield 

(R' ± ) 2 = -1, (15) 

which effectively changes the index v* by 2 and switches the topological charge between even 
and odd. Opposite parities naturally cancel out. 

As an example, consider a disclination loop with a reference profile of —1/2 at v = 0, 
that rotates by ip\ around its tangent, flips to a +1/2 profile around an axis at an angle 9\ 
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and with parity P\, then rotates by ip2 and flips back to —1/2 around an axis at 82 with 
parity P 2 . The total rotation of the profile along the entire loop becomes 



The exponent must be a multiple of 7r/2 to satisfy the continuity equation, which puts a 
restriction on the angles 0^2 and #1,2- This generalises the third-integer quantization of the 
self-linking number. When this condition is met, the above expression can be compared to 
k u * to extract the index v* . 

Note that the self-linking rotations of +1/2 and —1/2 parts have opposite-signed contri- 
butions to the index v* . Taking a different part of the loop as a reference cyclically permutes 
the order of the transformations. For example, if a +1/2 part serves as a reference, the en- 
tire exponent changes sign. This does not affect the R' = ±1 rotations {y* G {0,2}), but 
exchanges the signs of R' = ±k for linked disclinations. The same effect is observed if the 
parametrization of any of the disclination loops is reversed, which also reverses the linking 
numbers in equation ([7]). This ambiguity stems from the requirement of the classification 
with homotopy groups to have a fixed reference - a base point [B] . The choice of reference 
is very important once more than one loop is present in the system. 

Most of the disclination loops can be analysed directly by the described formalism. Wher- 
ever the normal of the circle traced by the director on the unit sphere is not perpendicular to 
the disclination tangent, it can be simply regarded as being a— l/2ora+l/2 profile and the 
director can be unambiguously combed to the perpendicular plane. The decomposition thus 
takes advantage of easily recognisable physical features of the disclination profiles and can 
be counted by hand if provided with a properly visualised director field. In singular cases, 
such as a twist disclination line, which is exactly half-way between +1/2 and —1/2 profiles, 
appropriate quaternion expressions can be quickly retrieved using a procedure similar to the 
one described above. 

V. NEMATIC SURFACE TEXTURES 

To present above statements from a different perspective, we must first shortly review the 
theory behind nematic point defects. The topological charge of a set of defects contained 



R' =(-P 1 P 2 )oexp(-ke 2 )exp(0<f )oexp(-k^) exp(k^) 




(16) 
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FIG. 3: The correspondence between subsurface disclination loops and geometric features of 
nematic textures on a sphere- like surface (a). The patches of perpendicular director (red for 
outgoing and blue for ingoing field) are segmented by boundaries (orange) with the director parallel 
to the surface. Boundary signature s = ±1 corresponds to +1/2 (b) or —1/2 (d) profiles of the 
subsurface disclination loops, while the grains - places where the director is parallel to the boundary 
- correspond to the twist profile, which is an intermediate state in transition between the other 
two elementary disclination profiles. 

inside a surface is determined by the second homotopy group class of the nematic texture 
on said surface. A convenient way of calculating the topological charge is introduced in 
reference [23]. It states that a surface texture can be split into patches (2-dimensional 
objects) where the director is perpendicular to the enclosing surface, separated by boundaries 
(1-dimensional objects) where the director is parallel to the surface, but perpendicular to 
the boundary itself, and grains (isolated points on the boundaries) where the director is also 
parallel to the boundary of which it is part. The topological charge of the texture can then 
be calculated by counting the contributions of boundaries t = ±1 and grains g = \. 

9=l + $> + X> (17) 
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The first term is a half of the Euler characteristic of the enclosing surface and is related to 
the fact that a radial hedgehog has a topological charge of q = 1. 

There is a striking similarity between the classification of disclination loops and the clas- 
sification of nematic textures. Both theories observe the behaviour of the director in a local 
frame, aligned to the surface or the curve in question. This splits the global contribution 
of the shape's topology from the local geometric contributions. For disclination loops, the 
global part is the Ro = — 1 factor that describes the nonzero topological charge of a simple 
Saturn ring defect, while for surface textures, it is the +1 contribution of the Euler charac- 
teristic that describes the nonzero topological charge of a radial hedgehog. Both approaches 
also assume the texture or the director profile can be smoothly transformed, so that all the 
changes happen in narrow transition zones. 

The Figure [3] shows that every surface texture can be alternatively interpreted as a 
homeotropic core with added disclination loops below the surface of interest at the places 
where the boundaries are, and with disclination profile transitions where the grains are. 
The condition that the grains come in pairs and change the topological charge by 1 if they 
have the same parity, corresponds exactly to the same rule for the profile transitions. The 
topological charge of the texture can therefore also be calculated by adding the core +1 
defect (the Euler characteristic) to the charges contributed by Janich's indices for all the 
virtual subsurface disclination loops below the boundaries, 



q 



l + n + ^gj mod 2, (18) 



for n boundaries and grains gj = ±|. One big difference between equations (17) and (18) 
is that the result is only given by modulo 2 in the latter. This is a natural consequence of 
dealing with line defects instead of point defects (see Ref. [IE]). However, the narrow subset 
of disclinations that can act as subsurface representations of boundaries, are restricted by 
the surface of the homeotropic core, so there must be no linking and self-linking. 

The correspondence between the textures and disclination loops is not necessarily only 
theoretical. If the surface on which we are observing the texture lies just above the disclina- 
tions, the geometric features on the surface texture will reflect the underlying disclination 
geometry. From this, another strong conclusion can be derived. The defect rank, a geometric 
measure for classifications of point charges, introduced in Ref. (23] , is only valid if there are 
no grains on the surface. By the correspondence principle, we can state that if there are no 
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disclination transitions, no linking and no self-linking in the system, the defect rank can be 
used and the usual intuition of assigning proper integer charges to both disclination loops 
and point defects is valid. As soon as the system becomes entangled or there are different 
disclination profiles present in the system, there are possible discrepancies in the full integer 
conservation of the topological charge and only the even and odd conservation law holds. 

VI. CONCLUSION 

The topological properties of nematics are an interesting topic, crucial for interpretation 
and design of experiments. Although the basic, strictly topological formalism has been 
known for a long time, easy application of the theory to the real-world examples requires 
geometrization of the basic topological concepts. By considering local director behavior, 
we can account for restrictions enforced by the energy minimisation and thus make a fine- 
grained classification using not only the homotopy theory but also additional geometric 
information. 

We have decomposed a general disclination loop into a sequence of local rotations that 
must add up continuously to a rotation that satisfies the director continuity in a closed 
disclination loop. Quaternions are an analytically efficient way of enumerating the local be- 
havior of the director profile. The use of quaternions further reinforces the close relationship 
between the uniaxial and biaxial nematics, as the unified handling of the entire cross section 
removes the differences between both phases. 

As an extension of the theory of —1/2 defects, the introduced formalism describes a 
wider range of disclination networks. Notable examples of systems that exhibit transitions 
between different director profiles are confined cholesteric phases [T^l [To] , optically induced 
defects [21] and transient defects that are created with quenching from the isotropic phase 
[25] . Studies of random blue phases [26] and entangled structures in random or structured 
pores [12] could also benefit from improved classification of disclination profile variations. 

The formalism we presented also shows a stunning connection with the texture decom- 
position, suggesting that the approaches of geometrisation for different kinds of defects can 
be unified, leading toward a consistent set of tools for enumeration of defects. 
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